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ABSTRACT 

We present environmental dependence of dusty star forming activity in and around the cluster 
RXJ1716. 4+6708 at z = 0.81 based on wide-field and multi-wavelength observations with 
the Prime Focus Camera on the Subaru Telescope (Suprime-Cam) and the Infrared Camera 
onboard the AKARI satellite (IRC). Our optical data shows that the optical colour distribution 
of galaxies starts to dramatically change from blue to red at the medium-density environment 
such as cluster outskirts, groups and filaments. By combining with the AKARI infrared data, 
we find that 15 ^m-detected galaxies tend to have optical colours between the red sequence 
and the blue cloud with a tail into the red sequence, consistent with being dusty star forming 
galaxies. 

TTt: ■ The spatial distribution of the 15 /zm-detected galaxies over ~ 200 arcmin 2 around the 

cluster reveals that few 15 /im galaxies are detected in the cluster central region. This is 
probably due to the low star forming activity in the cluster core. However, interestingly, the 
fraction of 15 /im-detected galaxies in the medium-density environments is as high as in the 
low-density field, despite the fact that the optical colours start to change in the medium-density 
environments. Furthermore, we find that 15 /.tm-detected galaxies which have optically red 
colours (candidates for dusty red galaxies) and galaxies with high specific star formation rates 
are also concentrated in the medium-density environment. These results imply that the star 
forming activity in galaxies in groups and filaments is enhanced due to some environmental 
effects specific to the medium-density environment (e.g. galaxy-galaxy interaction), and such 
a phenomenon is probably directly connected to the truncation of star forming activity in 
galaxies seen as the dramatic change in optical colours in such environment. 

Key words: galaxies: clusters: individual: RXJ1716. 4+6708 — galaxies: evolution — galax- 
ies: starburst — large-scale structure of Universe. 



> 

in 

o\ 

ON 

o 

oo 

o 

> 

X 



1 INTRODUCTION 

1.1 Galaxy properties as a function of environment 

Galaxies live in various environments. Recent redshift surveys have 
shown filamentary large scale structures in the local Universe. In 
the distant Universe, at least up to z < 1, similar filamentary nature 
of large scale structures is found around clusters through wide field 
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observations of distant clusters (e.g. lKodama et al.l2005h . There are 
also some hints that the large scale structure is present at much 
highe r redshifts up to z ~ 6 iShimasaku et alj[2003l ; lOuchi et"al] 
l2005t) . 

Environment must have played an important role in the history 
of galaxy evolution since galaxy prop erties are stron gly dependent 
on environment. This is first noted bv lDressleil J 198Ct) . who showed 
that early-type galaxies dominate in high density regions while 
late-type galaxies tend to live in low density regions. This trend is 
called "morphology-density relation". After the Dressler's work, 
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this inter esting trend was confirmed and extended by many au- 
thors (e.g.lPostman &Gellerl 19841 : iDressler et alj 19971 ; IGoto et al.l 



2003d: IPostman et al.l2005h . However, it is still unclear what is the 
key physical process to produce the morphology-density relation 
or other environmental dependence of galaxy properties. 

The morphology-density relation can be understood, at least 
partly, as a result of morphological transformation of a substantial 
number of galaxies when they entered high-density environments. 
Many mechanisms to suppress the star forming activity and to con- 
tribute to the morpho logical transformation have been proposed 
(see the review by iBoselli & Gavazzil l2006h . For example, ram- 
pressure stripping due to the interaction w ith hot plasma gas filled 
in the cluster core (e.g. lGunn & Gottl 19721) , is expected to be effec- 
tive in rich cluster cores. High-speed encounters between galaxies, 
which are often called "galaxy h arassment", should occur in very 
high-density environments (e.g. iMoore et alj 1 1996h . In addition, 
mergers or galaxy-galaxy i nteractions should also co ntribute to the 
galaxy transformation (e.g. iToomre & Toomre|[l972l) . Interactions 
with the cluster potential may also cause a tidal force when galax - 
ies pass the central region of clusters (e.g. lBvrd&Valtonerll99Ch . 
Another me chanism that can be effective is the so-called "strangu- 
lation" (e.g. lLarson. Tinslev, & CaldweHlll980D , which leads to a 
slow decline in star formation rate after a galaxy falls into a more 
massive (i.e. group or cluster) halo. Identification of the key pro- 
cesses behind the galaxy transformation is one of the major remain- 
ing issues in galaxy evolution. 

It is well known that the fraction of blue star forming galax- 
ies in clusters increases towards higher redshifts (Butcher-Oemler 
effect; iButcher & OemleJ 1 19841) . Therefore, the transition from 
blue active galaxies to red passive galaxies should be more com- 
monly seen in distant clusters. We can expect to see directly 
such trunc ation in action th rough observations of distant clusters 
(see also IGoto et alj|2003al) . However, importantly, it is reported 
that most of the actions take place in the outskirts of clusters 
rather than in the cluster core. In fact, some recent studies fo- 
cus on the galaxy properties in the surrounding regions of clus- 
ters and try to identify the environment wher e the truncation of star 
formation occurs in accreted galaxies ( e.g. I Abraham et al.|[l99^ ; 
iBalogh et alJl 19991 : [Pimbblet etai] |2002h . In such environment, it 
is reported that passive spirals (i .e. spiral morphol o gy but no star 
forma tion) tend to be found (e.g. lGoto et alj|2003bh . lKodama et alj 
d200ll) performed a wide-field imaging of the CL0939 cluster at 
z — 0.41 and discovered that the colour distribution changes 
dramatically at the intermediate density environment which corre- 
sponds to groups/fi laments. A very similar result was reported by 
iTanaka et alj d2005h for the surrounding regions of higher redshift 
clusters, CL0016 at z = 0.55 and RXJ0152 at z = 0.83. They 
suggest that the intermediate density environment such as groups 
or filaments around clusters are the very sites where the truncation 
of galaxies is taking place. These pioneer works are really telling 
us the need for wide-field observation of distant clusters in order 
to study the physical mechanisms that are responsible for the trun- 
cation of galaxies from active phase to passive phase during the 
course of hierarchical assembly of galaxies to clusters. 



1.2 Dusty star forming galaxies in the local and distant 
Universe 

It is well known that red galaxies tend to have little star forming 
activity (passive galaxies), while blue galaxies have on-going star- 
forming activity (star-forming galaxies) at a given redshift. How- 
ever, the classification of passive or star-forming galaxies based 



only on their optical colours is sometim es highly uncertain. In 
fact, lHaines. Gargiulo. & Merluzzil J200^) showed that ~ 30% of 
field red sequence galaxies selected from optical colour-magnitude 
diagrams have on-going star formation activity with EW(Ha) > 
2A, using t heir local galaxy s a mples from Sl o an Di gital Sky Sur- 
vey (SPSS: lYork et alj|2000h . iDavoodi et al] (12006)) also showed 
in their SDSS and SWIRE survey dLonsdale et aljT2003l) that ~ 
18% of their samples have red optical colours and infrared ex- 
cess at the same time, which include both AGNs and on-going 
dusty star-forming galaxies. Similar results are shown for cluster 
red sequence. For example, IWolf. Gray, & Meisenheimeil (2005) 
studied the Abell 901/902 cluster at z = 0.17, and showed 
that ~ 22% of red sequence galaxies are dusty red galaxies 
from the SEP fitting using the m edium-band survey (COMBO 17; 
IWolf, Meisenheimer, R6seill200ll) . 

Instead of using optical colours, optical emission lines (e.g. 
[Oll](A = 3727A) and Ha (A = 65 63 A)) are often and widely 
used to measure star formation rates (Kennicutt 1998). However, 
these lines, especially [Oil] lines in the rest-frame ultra-violet, are 
attenuated by inter-stellar dust in the galaxies. Also, [Oil] strength 
is affected by AGN contribution, if any, and dependent on metallic- 
ity as well. Although Ha line is a much better indicator than [Oil] 
in this respect (and is in fact one of the best indicators of star for- 
mation rates among emission lines), it is redshifted to near-infrared 
(NIR) regime at z Z 0.5, where large format wide-field camera or 
spectrograph has become available only recently. Given these dif- 
ficulties in optical and NIR observations, infrared (IR) luminosity 
of a galaxy that can be sensitively obtained by space telescopes, 
serves as an ideal measure of star format ion rate, and it is also 
well calibrated with local starburst galaxies J Kenni cutt 1998). Since 
the total IR luminosity of a galaxy can be estimat ed through sin- 
gle br oad-band imaging at mid-infrared (MIR) (e.g. lTakeuchi et all 
2005), it is critically important to observe clusters in MIR bands 
as well as in the optical/NIR in order to reveal the hidden star for- 
mation activity and hence trace its true star formation history. This 
is especially true at high redshifts, because the luminous infrared 
galaxies (LIRGs) which have lO 11 !,© < Lir(8 - lOOO^im) < 
lO 12 Z/0 and ultraluminous infrared galaxies (ULIRGs) which have 
£ir(8 — 1000/xm) > 10 12 Lq are more comm only seen in the dis- 
tant Universe than in the l ocal Universe (e.g. ISanders & Mirabel 
ll996l : lLe Floc'h et alj2005h . 



1.3 Infrared observation of galaxy clusters 

Taking these situations into account, wide-field MIR study of 
galaxy clusters covering entire structure around the cluster is re- 
quired to investigate the "true" environmental dependence of star 
formation activities of galaxies. However, until recently, infrared 
studies o f galaxy clusters ha ve been conducted mainly with the ISO 
satellite dKessleretal.1 19961) . which h ave been limited to very inner 
regio ns of clusters (see the review by Metcalfe, Fadd a. & Bivianol 
120051) . Even though, some of these studies showed the importance 
of large amount of hidden star f ormation activity in the cluster en- 
viron ment (e.g.[P uc et al. 2002). The recent advent of the Spitzer 
MIPS (Rieke et al. 2004]) has en abled us to investi gate wider fields 
of distant clusters up to z ~ 0.8. lGeach et al. (2006) observed very 
wide field (25' x 25') around CL0024 (z = 0.3 9) and MS0451 
(z = 0.55 ) clusters by m osaicing MIPS images. iMarcillac et alj 
1 20071) and lBai et alj J2007I) observed RXI0152 and MS 1054 (both 
at z = 0.83), respectiv ely, but their fiel d coverage is limited only to 
cluster central regions. iBai et al.N2007}) actually imply that for rich 
clusters at z ~ 0.8 even wider-field infrared studies are needed. 
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In the local Universe, it is well established that star for- 
mation activity depends strongly on environment in the sense 
that it systematical ly declines towards higher density regions (e.g. 
iGomez et alJl2003D . Recently, however, surprising results are re- 
ported where such a relationship between star format ion activity 
and l o cal galaxy density b ecomes inverted at z ~ 1 ( Elb az et"al] 
120071 : ICooper et al.l 120071) . i.e., star formation rate increases to- 
wards higher density environment at 2 ~ 1. This is naively ex- 
pected that one approaches to the formation epoch of cluster galax- 
ies which is probably systematically skewed to higher redshifts 
compared to the formation epoch of field galaxies. Therefore, look- 
ing back the environmental dependence of star formation activity in 
galaxies as a function of redshift is a basic but vital step towards un- 
derstanding the environmentally dependent galaxy formation and 
evolution. 

In this paper, we conduct for the first time a panoramic MIR 
study of a z ~ 0.8 cluster, over a very large area including sur- 
rounding groups, filaments and outer fields. Throughout this paper, 
we use Om = 0.3, = 0.7, and H = 70 km s _1 Mpc _1 . Mag- 
nitudes are all given in the AB system, unless otherwise stated. 
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2 DATA 

2.1 RXJ1716 cluster 

The RXJ17 16.4+6708 cluster at z = 0.81 that we study in this 
paper was first discovered in the ROSAT North Ecliptic Pole Sur- 
vey (NEPilHenrv et al]l997t) . Optical spectroscopy was performed 
bv lGioiaetaiT Tl999) and 37 cluster members were identified. Us- 
ing these spectroscopic samples. iGioia et aT. I GUI) determined the 
velocity dispersion of this cluster to be a = 15221^q km s _1 . 
This value is relatively large for its rest-frame X-ray luminosity of 
L bol = 13.86 ± 1.04 x 10 44 erg s _1 and the t emperature kT = 
6.8lj g keV which are ba s ed on Chandra data tettori et al.ll2004l 
see a lso IGioia et al.l 1 19991 IVikhlinin et al.l 120021 and iTozzi etall 
2003). Therefore, it is suggested that this cluster is not totally viri- 
alised yet. In fact, the brightest cluster galaxy (BCG) o f this cluster 
is loc ated on the northwestern edge of the structure dClowe et al.l 
Il998h . which may be linked to the fact that the structure of this 
cluster is still being formed. The weak-lensing mass of this cluste r 
is estimated to be 2.6 ± 0.9 x 10 14 /!, _1 M Q dClowe et alJI 19981) . 
This is consistent with the estimated mass based on the the X-ray 
data in lEttori etall fe004l) ; M tot = 4.35 ± 0.83 x 10 14 M Q . As 
the X-ray properties of this cluster suggest, this cluster is relatively 
rich for a cluster at this redshift. 

It has been known that this cluster has a small subcluster or 
group to the northeast of the main cluster, and the morphology of 
the X-ray image of this c luster elongates towa rds the subcluster 
(e.g. |jeltemaet al1l2005l) . iKovama et all (2007) performed wide- 
field optical imaging of this cluster and discovered prominent large- 
scale structures penetrating the cluster core and the second group 
of this cluster towards the southwest of the cluster core based on 
the photometric redshift technique (see Section l2T2t . 



2.2 Optical data 

We h ave been conducting the PISCES project (Koda ma et al.l 
120051) . Taking advantage of the wide field coverage of the 
Prime Focus Camer a on the Subaru Telescope (Suprime-Cam; 
Miva zaki et"afl f2002), we have investigated various environments 
around the X-ray detected distant clusters at 0.4 <: z Si 



Figure 1. The response function of the AKARI L15 filter and the template 
SED of starburst galaxies with 10 Lq from Lagache et al. (2004). Note 
that our L15 filter neatly captures the PAH broad emissions (7.7 (im and 
8.6 fim in the rest-frame) from z ~ 0.8 galaxies. 



1.3 (e.g. [Kodama et alJl2005l: iTatiaka et alJl2005L 2006, 2007a,b; 
iNakata et alJl2005l:lKovama et alj|2007ir As a part of this project, 
IKovama et al.l ( 120071) conducted a deep and wide-field optical study 
of the RXJ1716 cluster with the Suprime-Cam. We observed this 
cluster in Vffi'z'-bands which neatly bracket the 4000A break fea- 
ture of z ~ 0.8 galaxies. Catalogues were created with z' < 24.9 
galaxies, which corresponds to 5a detection limits of our z'-band 
image. We select the cluster mem bers by photometric redsh i ft tech - 
nique (phot-z) using the code of IKodama. Bell, & Bowel d 19991) . 
and showed the prom inent large-scale stru ctures around the cluster 
(see Figs. 3 and 4 of IKovama et alj|2007l) . Galaxies with 0.76 < 
Zphot < 0.83 were used to map out the large-s cale structures 
and th is member selection is used in this paper. See IKovama et al.l 
(2007) for the summary of our optical data, data reduction, and the 
combined colour image of this cluster. 



2.3 Infrared data 

We obtained deep and wide-field IR imagin g data for the RXJ1716 
cluster with the Infrar ed Camera (IRC: | Onaka et alj|2007l) onboard 
the AKARI satellite (Mur akami et alj|2007t ). Due to the good vis- 
ibility of NEP directions from AKARI, we could get very deep 
data for the RXJ1716 cluster which is located near the NEP (see 
Section |2~TV The observations were executed from November 26 
to December 8 in 2006, on an open-use programme CLNEP(PI: 
T. Kodama) and a mission programme CLEVL(PI: H.M. Lee). 
We observed this cluster in the N3(3.3 /xm), S7(7.0 /xm) and 
L15(15.0 fim) filters, using the Astronomical Observation Tem- 
plates (AOT) for deep observation (IRC05) with an AOT filter com- 
bination parameter of "b" (see AKARI IRC Data User Manual, 
lLorente et alj|2007l) . We should note that the L15 fi lter neatly cap- 
tures the polycyclic aromatic hydro carbon (PAH; IPuget & Leged 
119891) emissions from star-forming galaxies at z ~ 0.8 (Fig.[T). 
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Figure 2. The three AKARI FoVs (Fl, F12 and F2) on the z'-band image (lej 
in both of the figures. The size of each panel is ~ 18' X 18'. 

We set the three field of views (Fl, F12 and F2) to cover the 
large scale str uctures around the RXJ1716 cluster photometrically 
discovered bv lKovama et alJfeOOTT) (see Figs.[2]and[7p). The num- 
ber of the pointed observations and corresponding integration time 
for each field are summarised in Table[T] Note that 10 and 30 frames 
with effective exposure times 44.4 and 16.4 seconds were included 
in a single pointed observation for N3 and S7/L15, respectively. 
Therefore, typically, one pointed observation corresponds to 7.4 
min and 8.2 min exposure observation for N3 and S7/L15, respec- 
tively. IRC has a wide field of view (about 10' x 10'), and our 
observed fields Fl, F12 and F2 covered about 200 arcmin 2 in total 
(see Figs.[2]and|7jj). We stress here that our study is the first attempt 
to cover such a wide field around the z Z 0.8 cluster in MIR bands. 
We can study the cluster core to outskirt regions at once and this is 
critically important to discuss the environmental dependence of the 
infrared properties of galaxies in and around high-z clusters. 

The data w ere reduced with th e IRC imaging pipeline (ver- 
sion 20071017. |OnakaetalJ 120071) in a standard manner (see 
lLorente et al.ll2007l) . We ran the pipeline for each of the pointed 
observational data, and we made the reduced image for each of 
the pointed observations at first. In this process, we subtracted the 
median filtered sky with 20 x 20 kernel for individual frame be- 
fore co-adding the frames. Note that the pipeline often fails to 
align L15 frames because of the small number of firmly detected 
sources. Therefore, we used "coaddLusingS" option in the pipeline, 
which us es the alignment of the S7 frames to align the L15 frames 
(see also lLorente et alj|2007l and lWada et~al]|2007h . After that, we 
matched the position of all the images generated from each pointed 
observational data using the IRAF tasks 'geomap 7 and 'geotran' , 
and we co-added all the images together. The final 15 /im image is 
shown in the right panel of Fig. [2] 

Since the exposure times are not uniform over the observed 
field, limiting flux depends on the position. However, as we ob- 
served with 14 pointings for Fl field in N3 and S7 band and with 
15 pointings for F12 field in L15 band (see Table[]}, the depths of 




and final co-added L15-band image (right). North is up and east is to the left 



the data depend mainly on whether a position is covered by these 
deep observations or not. We thus divide each image into "deep" 
and "shallow" regions and estimate the sky noise limits of both 
regions by a simple aperture photometry of randomly distributed 
apertures in all the images. Note that we do not use areas near 
very bright sources. The aperture diameters were set to 8", 10" 
and 11" for N3, S7 and L15, respectively. These correspond to 
twice the size of FWHM of the point spread function of each im- 
age. We also used these diameters for source extraction (see be- 
low). The 5a limiting fluxes of deep and shallow regions for each 
band is su mmarised in Table [2] Sources were extracted using SEX- 
TRACTOR jBertin & ArnoutslTl996l) . Using the same size of aper- 
tures (i.e. twice the size of FWHM of each image), we measured 
FLUX_APER of sources in each image. We made the 3 ^m catalogue 
with FLUX_APER > 7.0/xJy in N3 band, the 7 fj,m catalogue with 
FLUX.APER > 21.9/LtJy in S7 band and the 15 catalogue with 
FLUX_APER > 66.5/xJy in L15 band. These 5a detection criteria 
for the deep regions approximately correspond to ~ 3.5cr detec- 
tion for the shallow regions (Table O. Also, we use FLUX_AUTO 
and MAG_AUTO as the total fluxes and the total magnitudes of the 
sources as shown in IWada et al.l < [2007l) . We note that all the scien- 
tific quantities such as total IR luminosities and star formation rates 
are measurements based on these total fluxes (see also Section [3~2l . 



3 ANALYSIS 

3.1 Cross identification between optical and IR sources 

We need to cross match between the AKARI IR sources and the 
Subaru optical sources. We should be very careful when we iden- 
tify the sources in a crowded region like a galaxy cluster. We have 
cross-identified the 15 /im sources in the optical (z'-band) image 
and created "15 fim member" catalogues in the following way. 

Firstly, we search for any optical counterpart(s) using a 8" ra- 
dius from each 15 fim source. This radius is sufficiently larger than 
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filter 


Fl 


F12 


F2 


N3 


14 (104 mm) 


2(15 mm) 


3 (22 min) 


S7 


14(115 mm) 


3 (25 mm) 


3 (25 min) 


L15 


3 (25 min) 


15 (123 mm) 


4 (33 mm) 



Table 1. Summary of the AKARI(IR) observation. The number of point- 
ings and corresponding exposure time for each filter and each field is shown. 



filter 


deep (5cr) 


shallow (5cr) 


N3 


7.0 /Jy 


11.9/iJy 


S7 


21.9 ixSy 


33.0 /iJy 


L15 


66.5 /tjy 


96.5 /jjy 



Table 2. 5a limiting fluxes of N3, S7 and L15 images. Deep regions ap- 
proximately correspond to Fl field of N3/S7 images and F12 field of L15 
image. Shallow regions are defined as the rest of the deep regions (see text). 

the FWHM of the PSF in the 15 fim image (~ 5"5). The relative 
positional accuracy between the z'-band image and the L15 image 
is also sufficiently smaller ( < 1") than this search radius. We find 
that 149 of 15 fim sources have at least one phot- z member galaxy 
(i.e. with 0.76 < z p i lot < 0.83 as used in iKovama et alj (2007) 
to trace the large scale structure) within this radius. Secondly, we 
carefully examine all of them by eye. We basically select the near- 
est object as the optical counterpart of each 15 fim source. If the 
optical counterpart is confirmed to be an isolated phot-z member 
galaxy, we put it in the "resolved 15 fim member catalogue" (see 
an example in Fig. [3^). On the other hand, due to the large FWHM 
of the PSF of the 15 /im image compared to the optical image, some 
15 fim sources are not resolved in the source extraction process al- 
though they show multiple peaks or extended shape in the 15 fim 
image. We find that the number of such extended sources are not 
negligible. Therefore, in such a case, we carefully look at the shape 
of the 15 /im image by eye and judge if the 15 fim multi-peak or 
extended source is associated with any phot-z members. If we find 
any phot-z members associated with the 15 fim source, we put it 
into the second catalogue, named "unresolved 15 /im member cat- 
alogue" (see the examples in Fig.[3p,c). 

Most of the "resolved" 15 /im members have optical counter- 
parts within 2" radius (median value is 0"7), while counterparts 
of "unresolved" 15 /im members are found at relatively larger dis- 
tance (median value is 2"5). We note that in the latter catalogue 
there are some 15 /im sources which have more than one optical 
counterparts (Fig.[3p,c). For such sources, there is a possibility that 
the 15 /im flux of them is measured as a sum of the fluxes from 
more than one objects. Therefore, we cannot know its separated in- 
dividual 15 /im fluxes correctly and we can only know an upper 
limit. The cross identification itself is more difficult for unresolved 
15 /im members due to the complexity of the shape of the 15 /im 
sources. Given these problems, we have to treat them carefully in 




Figure 3. Examples of the 15 /im member galaxies. The panel (a) shows 
an example of "resolved" cluster member while the panels (b) and (c) show 
"unresolved" cluster members. In each panel, the large circle indicates the 
position of the detection of 15 fim sources, and its radius is 8". Thin-line 
small squares indicate the positions of the optical sources, and the size of 
them is 3" X 3". Thick-line squares show the positions of phot-z selected 
galaxies. In the case of the (b) and (c) panels, we select three and one phot-z 
members as counterparts of the 15 fim source, respectively. 

the following sections, and in fact we do not include these "unre- 
solved" 15 /im members in the following analyses whenever the 
15 fim fluxes are required. We count these sources only if detection 
at 15 fim is concerned. However, we have confirmed that the inclu- 
sion of such "unresolved" sources or not does not affect our results 
(see also Section 4.4). 

For the total of 436 sources detected at 15 fim in our AKARI 
field, we find that 91 of 15 fim sources are associated with phot-z 
members. In this sample, 54 sources are turned out to be resolved 
15 fim members. The remaining 37 sources are unresolved cases 
and we have 65 unresolved 15 fim members. This is because about 
a half of these 37 sources are associated with more than one phot-z 
members. Therefore, we have 54 resolved and 65 unresolved 15 fim 
members in total as the numbers of optical counterparts. For the re- 
solved 15 fim members, we search for a counterpart in N3 and S7 
images in the same way as the cross identification between L15 
sources and optical sources described above (i.e. after searching 
for N3/S7 counterparts using 8" radius from each 15 fim member, 
checking them by eye ). In the resolved 15 fim members, 18 sources 
are identified in both N3 and S7 bands, and 24 sources are identi- 
fied only in N3 band. The rest of them are not well distinguished as 
a single object in either N3 or S7 bands. Our main interest is in the 
15 fim detected cluster members because these sources are the can- 
didates for dusty star forming galaxies. Therefore, we use 15 fim 
detected cluster members in the following sections regardless of the 
identification in N3 or S7 images. 

Note that photometric redshift is not always the precise crite- 
rion to select the cluster member galaxies. There remains a possi- 
bility that some fore-/background contaminations are included in 
our 15 fim member catalogues. Unfortunately, we cannot estimate 
the probability that a galaxy out of the 0.76 < z p hot < 0.83 range 
falls into this redshift range by our phot-z estimation because we 
do not have spectroscopic information of normal galaxies in our 
RXI1716 field. However, we do know some spect roscopic r edshifts 
for the X-ray sources in our RXJ1716 field from lKim et alj J2007t) . 
We find that only 1 out of 27 such sources is miss-identified as 
our phot-z member. Therefore, the fore-/background contamination 
in our phot-z samples is expected to be very small, although firm 
conclusion waits for a similar check using normal galaxies, which 
cannot b e done at the moment because of unavailability of their 
redshifts. iMarcillacetailfeOOTl) , who studied the RXJ0152 cluster 
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at z — 0.83, estimated that the fraction of 151 spectroscopically 
confirmed non-members that fall into the 0.76 < z v hot < 0.88 
ra nge is only ~ 7%. T hey used the photo metric redshifts estimated 



Tanakaetafl l l2005l) e. phot-z code in Kodama, Bel l. & Bowed 
using the Suba ru VRi'z' data), wh ich is the same as that 
used in this study. iMarcillac et ai] (2007) also showed that the 
photometric redshifts worked well for their spectroscopically con- 
firmed MIR member galaxies (\Az\ = \z p hat — z spec \ ~ 0.01). 
Therefore, the photometric redshift selection would not do any crit- 
ical harm on the statistical properties of our 15 fim members, ei- 
ther. We should, however, confirm the physical association of these 
sources to the cluster through spectroscopic follow-up observations 
in our future work. 



3.2 Derivation of star formation rates for the 15 /im 
members 

The AKARI's L15 band corresponds to ~ 7 — 9 /im in the rest 
frame of the RXJ1716 cluster at z — 0.81. This wavelength range 
neatly includes 7.7 and 8.6 broad line emission features of 
PAHs (Fig. QJ. It is well known that MIR broad band luminosity 
correlates with total infrared luminosity (Lis.) and hence star for- 
mation rate t hrough a good correlation between Lir and star for- 
mation rate fcemicutd[l998t) . It is particularly true in this case, 
since the PAH emissions originate from photo dissociation regions 
associated to star forming regions and the intensity of PAH fea- 
tures thejns^ly£s_arejn_good correlation with star formation rate 
(e.g. lCharv&Elbazll200lh . 

To derive Lir of each galaxy, we use the template SEDs of 
starburst galaxies from Lag ache et al.ld2004h . We calculate the con- 
version factor of our 15 fim (rest frame ~ 8 fim) luminosity to 
the total (8-1000 /im) luminosity using their 5 template SEDs (i.e. 
templates for 10 9 , 10 10 , 10 11 , 10 12 and 10 13 L Q starburst galax- 
ies). The derived correlation between vLv (8 /im) (AKARI L15 
band flux) and Lir (8-1000 ^m) is shown in Fig. [4] Using this 
correlation, we estimate the total IR luminosity from the measured 
total 15 fim flux (i.e. FLUX.AUTO value from SEXTRACTOR ) of 
each 15 fim member. Then, using the Kennicutt (1998) calibration: 



SFR(M Q • yr" 1 ) = 4.5 x 10" 



44 T I -In 

Lir (erg • s ) 



we derive its star formation rate. A histogram of Lir and corre- 
sponding SFR in our total samples are shown in Fig. [5] Due to the 
flux limit of our MIR observation, almost all of our MIR samples 
are turned out to be LIRGs with Lir > 10 1 Lq. Our 15 fim flux 
limit of 67 fiSy corresponds to SFR of ~ 25 Mo/yr. We calculate 
Lir and SFR only for the "resolved" 15 fim members since we 
cannot separate and measure individual 15 ^m fluxes of the unre- 
solved sources (see Section l3Tt . We should note that we do not 
require accurate absolute values of Lir and SFR in this paper and 
our conclusion will rely rather on their relative values as a function 
of environment. 

There is an issue that our IR galaxy samples may contain con- 
tamination from active galactic nuclei (AGNs). It is reported that 
the mean AGN fraction in clusters is only ~ 1% from optical spec- 
troscopic surveys, (e.g. Dressi er et al ] | 19991 for their 10 clusters at 
0.37 < z < 0.56). However, an exce ss of X-ray detecte d AGNs 
is also reported in some clusters (e.g. iMartini et alj|2002l) . In our 
15 fim member galaxies, we identified three X-ray point sources 
using the Chand ra X-ray point source catalogue constructed by 
iKim et al.l d2007h . Although we admit that we cannot detect all 
AGNs in X-ray, these three sources are strong candidates for AGNs 




log (yhv, 8 /urn) (L^ 



Figure 4. The relation between uLu(8fim) from AKARI data and Lir (8 — 
1000/im) derived from the template SEDs of Lagache et al. (2004). Five 
points correspond to their five SED templates of starburst galaxies. The 
AKARI 15 (im flux from a z = 0.81 galaxy can be converted into a total 
IR luminosity using this relation (see text). 
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Figure 5. Distribution of infrared luminosity Lir (8-1000 fim) and corre- 
sponding SFR for the 15 fim cluster members. Note that Lir and SFRs are 
calculated only for the resolved 15 fim members. 
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Figure 6. (a): The colour-magnitude diagram in R — z 1 vs. z' for all the member galaxies in the field covered by AKARI FoVs. Open circles and triangles 
represent the resolved and the unresolved 15 (im members, respectively. The three sizes of the symbols (large, middle and small) indicate 15 fim flux with 
(/(15/im) > 160/iJy, /(15/im) > 100/iJy and /(15/im) > 67/iJy), respectively. The solid line shows the best-fitted colour-magnitude relation defined 
in Koyama et al. (2007), and the two long dashed lines indicate ±0.2mag from the best-fit colour-magnitude relation. The vertical and slanted short dashed 
lines show the limiting magnitudes in z'-band and ij-band, respectively, (b): R — z' colour distribution of all the member galaxies (open histogram) and the 
15 fim members (hatched histogram). Histograms are constructed using only the optically bright z' < 22.5 galaxies in the field covered by AKARI FoVs. 
The solid-line locus indicates the fraction of the 15 fim members in each colour bin (see the label and the tick marks on the right side). Error-bars indicate the 
Poisson errors. The open histogram shows two peaks which correspond to red/blue bimodal populations, while the hatched histogram shows a single broad 
peak located between the two peaks of the open histogram. 



and we consider that at least a part of their 15 fim fluxes are emit- 
ted from AGNs. Another possible technique to distinguish AGNs 
from starburst galaxies is to use rest-frame NIR colours, based on 
the different NIR SED properties of these two populations (e.g. 
IWebb et aiT 2006). Rest-frame NIR SEDs for starburst galaxies are 
relatively flat, while AGNs produces power-law SEDs which in- 
creases towards longer wavelength. In our data set, the slope of the 
NIR SEDs calculated using the fluxes in N3 (rest-frame 1.8 /im) 
and S7 (rest-frame 3.9 fim) would be useful. We therefore investi- 
gated the N3— S7 colour of each resolved 15 fim member galaxy to 
see if there is any possible candidate of AGN. In our 54 resolved 
15 /im members with /(15/im) > 67/iJy, 18 galaxies are identi- 
fied in both N3 and S7 bands, 24 are identified only in N3-band, 
and the remaining 12 are identified in neither N3 nor S7 bands (see 
also Section [3~TV Out of the 18 galaxies detected in all bands, only 
two have NIR colour that is redder than N3— S7=0.0, which could 
be considered as AGN candidates. However, we notice that these 
two sources are not detected in X-ray. 

In these ways, we have identified five AGN candidates in to- 
tal from our 15 fim member galaxies. However, we do not exclude 
these objects from our sample because we cannot know if all the 
MIR emissions from these sources are from AGNs. It is still possi- 
ble that starburst is also in place in the AGN host galaxies. We stress 
here that the effect of including these AGN candidates for our con- 
clusion is negligible (see more detailed discussion in Section |4~4l . 
We should note, however, that we do not have any information of 
NIR SEDs for the objects that are not detected in either N3 or S7 
bands. There remains some possibility that some of these objects 
could be associated with AGNs, but we cannot exclude them under 
present conditions. In any case, future spectroscopic follow-up ob- 



servations of the 15 /im members is necessary to discuss the AGN 
contamination further. 



4 RESULTS 

4.1 15 fim members on the colour-magnitude diagram 

In this subsection, we present optical properties of the 15 fim de- 
tected galaxies in the RXJ1716 cluster, mainly focusing on optical 
colours of these galaxies. In Fig.|6ja), we show a colour-magnitude 
diagram in R — z vs. z 1 for the entire field covered by the AKARI 
observations. All the cluster member galaxies selected on the basis 
of photometric redshifts are shown. The 15 fim cluster members 
are represented by open symbols. The "resolved" and "unresolved" 
15 fim members are indicated by open circles and triangles, respec- 
tively. Most of the 15 fim members are optically bright galaxies 
with z' g 22.5 mainly because of the limited depth of the MIR 
observations. 

The 15 fim members are distributed in slightly bluer side of 
the red sequence, which is often called a "green valley" region. We 
show in Fig.|6jb) the histograms of R — z' colour distribution for 
all the phot-2 selected cluster member galaxies with z' < 22.5 
(open histogram) and for only the 15 fim member galaxies with 
z < 22.5 (hatched histogram). A clear bimodality is seen in the 
open histogram, while the hatched histogram shows a unimodal 
distribution that fills the gap between the two peaks of red and 
blue galaxies. This trend is quantified by calculating the fraction 
of MIR detected galaxies as a function of 7? — z' colour using the 
z < 22.5 galaxies (a locus in Fig. [6fe). A similar colour distribu- 
tion was reported in lGeach et al.l ( I2006h for two z ~ 0.5 clusters 
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Figure 7. (a): R — z' colours of phot-2 member galaxies (0.76 < z(phot) < 0.83) as a function of local projected number densities of member galaxies 
within the area covered by our AKARI observations. The vertical dotted lines at log = 1.65 and logS 5i h = 2.15 define the "low", "medium" and 

"high" density environments as indicated. Three dashed lines represent the loci of the 25th, 50th and 75th percentile colours. The solid line shows the fraction 
of red galaxies in each environmental bin. Error bars represent the Poisson errors, (b): Spatial distribution of the phot-z selected galaxies in high/medium/low 
density environments. Large filled circles, open squares and small dots represent galaxies in the high, medium and low density environment, respectively. The 
contours show the local 2D number density of galaxies at 2, 3, 4, 5cr above the peak of the density distribution. We apply Gaussian smoothing (sigma= 0.2 
Mpc in physical) on each galaxy and combine the tails of Gaussian wings to measure the local density at a given point. A bin size of 0. 1 Mpc (physical) is used 
to draw iso-density contours. The three large dashed-line squares show the three FoVs of our AKARI pointing observations. Note that the medium-density 
environment corresponds to the outskirts of the cluster core, subclumps and filaments. 



based on Spitzer MIPS observations. The medium colours of the 
MIR detected galaxies bridging the "green valley" would be due to 
dust reddening of blue star forming galaxies. We also find that this 
trend can be seen even if we limit the sample for relatively bright 
(z' < 21.5) or faint (21.5 < z' < 22.5) galaxies. We should note, 
however, that this trend might be produced partly due to the small 
systematic difference in luminosity of galaxies between those in the 
green valley (brighter) and those in the blue cloud (fainter). 

It is interesting to note that some red-sequence galaxies are 
detected at 15 /im. It can be interpreted that they are forming stars 
but are highly reddened by dust. Recently, it is reported that the 
fraction of such MIR detected star forming galaxies on the red se- 
quence i n the total star forming galaxies i ncreases towards distant 
clusters iSainto nge. Tran. & Holden 2008). We hereafter call these 
optically red 15 members "dusty red galaxies", and discuss 
them further in Section|4~4l 



4.2 Optical colour transition in the outskirts of the cluster 

In this subsection, we define galaxy environment on the basis 
of local projected number density of galaxies. We use a nearest- 
neighbour density, which is widely used by many authors (e.g. 
Tanaka et al. 2005). The number density at each galaxy's position 
is calculated using the neighbouring galaxies located within the cir- 
cle of the radius that is equal to the distance from the galaxy to the 
5th-nearest galaxy (hereafter, Estjj). In this calculation, all the op- 
tically selected phot-z members are used. 

Fig- H a ) shows the R — z colours of individual member 
galaxies as a function of the above defined local density (E^th)- 



The colour distribution changes dramatically at log Esth ~ 2.0 
as the median (50 %) colour locus suggests. It can also clearly be 
seen in the fraction of red galaxies that starts to increase sharply 
at logEsth ^ 1.7. We divide the galaxies into three environmen- 
tal bins, namely, low-density (logSsth < 1.65), medium-density 
(1.65 < log E 5th < 2.15) and high-density (log T, 5th > 2.15) re- 
gions. The medium-density region is defined to a relatively narrow 
range of local density where the optical colour distribution starts to 
change. These definitions of the environments may seem a little ar- 
bitrary, but our results do not change if we slightly change the def- 
initions of the environments. In Fig.|7Jb), we show the spatial dis- 
tribution of the member galaxies in each environment with differ- 
ent symbols. The high-density region corresponds approximately 
to the cluster core region. The medium-density environment trace 
outskirts of the cluster core, groups and filaments. The low-density 
environment corresponds to the remaining fields outside the groups 
and filaments. We note that the median value of log E5 for all the 
galaxies with 0.76 < z(phot) < 0.83 in our entire Suprime-Cam 
field (34'x27') is ~ 1.3. Excluding the cluster region would only 
lower this value by a negligible amount. Therefore, we can reason- 
ably consider that our low-density region corresponds to general 
field at the z ~ 0.8 Universe at least in terms of the local density 
of galaxies. 

The above finding that the colour distribution dramatically 
changes in the "me dium-density" e nviron ment is consistent with 
the previous works. iKodama et ah 1 2001 ) show ed similar results 
for the z ~ 0.4 cluster. and lTanaka et al. (2005) also found simi- 
lar trends for z = 0.55 and z — 0.83 clusters, both based on the 
panoramic imaging with Suprime-Cam on the Subaru Telescope. 
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Figure 8. (a) Spatial distribution of the 15 /im members around the cluster for the entire field covered by our AKARI observation. Small dots show all the 
phot-z member galaxies. Open circles and triangles indicate the resolved and the unresolved 15 fim members, respectively. The sizes of the open circles and 
triangles indicate 15 ftm fluxes as in the previous figure (Fig. [6} . Contours are drawn in the same way as in Fig.|7p. The three tilted dashed-line squares show 
our three AKARI FoVs. (b): Close-up view of the central 5' X 5' region of the left panel. The meanings of the symbols are the same as in the left panel, 
but contour levels are changed to 3, 5, 10, 15, 20<r above the peak of the density distribution. In this plot, 15 ftm members are shown down to 3.5<r level as 
thin-line tiny open symbols. Note that few 15 fim members are detected in the very central region (i.e. <L 1' from the cluster centre). 



The current result for another cluster (RXJ1716) lends further sup- 
port to the scenario that galaxies start to be truncated in groups or 
outskirts of clusters before they entered a very high-density envi- 
ronment such as cluster cores. This scenario requires some mech- 
anism that works efficiently in relatively low-density environment 
and truncates star formation activity. 



4.3 Low star formation rate in the cluster core 

We show in Fig. [8ja) the spatial distribution of 15 (im member 
galaxies on top of the distribution of phot-z selected member galax- 
ies over the 18' x 18' area centred on the RXJ1716 cluster. This area 
corresponds to ~ 15 Mpc x 15 Mpc in comoving scale at z = 0.81. 
The open circles and triangles indicate the positions of the resolved 
and the unresolved 15 ftm members, respectively. The 15 fim fluxes 
of the 15 fim members are indicated by three different sizes of the 
symbols (brighter objects are shown in larger symbols). 

A close-up map of the central 5' x 5' region is shown in 
Fig. Ob). It is worth noting that a 5' x 5' field corresponds to the 
field of view of a single pointing with Spitzer MIPS, and our total 
spatial coverage with AKARI is ~ 8 times wider than that. Our 
data is deepest at this cluster central region because this region is 
covered by all the Fl, F2 and F12 fields. Therefore, we here plot the 
15 fim members all the way to the faintest end at 15 fim (> 3.5a de- 
tection). However, we rarely see 15 fim members in the very central 
region of the cluster (i.e. < 1' from the centre), although the num- 
ber density of the phot-z cluster members is the highest in this very 
central region. In contrast, the outer regions just outside this zone 
of avoidance of the 15 fim members, many 15 fim members are de- 
tected. In particular, they are distributed preferentially at the region 
~ 1' north-east and south-west away from the centre. These regions 



are located right on the large scale NE-SW filament penetrating the 
cluster core and hosting two distinct groups (see Fig. [8^). 

To show it more quantitatively, we plot the cumulative radial 
profiles of the member galaxies as a function of distance from the 
cluster centre out to 4' in Fig. [9] Fig.|9]clearly shows that the 15 fim 
members are less concentrated in the < 1' region of the cluster than 
the general cluster member galaxies. From the slope of the curve for 
the 15 fim members, we see that many 15 fim members are located 
in the radius range of ~ 1.5-2.0 arcmin from the centre, which in 
fact corresponds to the region just around the cluster centre (see 
also Fig.(8j>). A Kolmogorov-Smirnov test on these two subsamples 
shows that the probability that these two are from the same parent 
is only ~ 3.5 %. A more detailed analysis of the distribution of the 
15 fim member galaxies will be carried out in Section |4"l4l 

4.4 Enhancement of dusty star formation activity in the 
cluster outskirts 

For each environment (i.e. low, medium and high density), we cal- 
culate the fraction of the 15 fim members in all the member galax- 
ies with z < 22.5, including both resolved and unresolved 15 fim 
members. We show in Fig.[l0]the fractions of the 15 fim members 
with /(15ftm) > 67fJy for red, blue and red+blue cluster member 
galaxies as a function of the environment. The boundary between 
red and blue is set as 

{R-z) = 2.51 - 0.049 x z, 

w hich is defined as 0.2 mag bluer than the best-fitted CMR shown 
in lKovama et alj j2007t) (see also Fig.|6^). We only used the galax- 
ies with z 1 < 22.5 because the number of 15 fim members with 
z' > 22.5 is very small due to the depth of the 15 fim image (see 
Fig.ISi). 
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We show in Fig. [10] that the 15 fim fraction for all galaxies 
decreases dramatically in the high-density environment compared 
with low and medium-density environments. This represents the 
very low star forming activity in the cluster core. On the other hand, 
we find that the 15 /im fraction is still high in the medium-density 
environment compared with the low-density fields, although we 
showed in Section 14.21 that the galaxy colours starts to sharply 
change in the medium-density environments. Although the statis- 
tics is not very good, the fraction seems to be higher in the medium- 
density environment than in the low and high-density environ- 
ments. Such a trend is likely to be strong for the red galaxies. This 
may indicate that in the medium-density environments dusty red 
galaxies are preferentially produced and/or the star forming activ- 
ity is once enhanced for some galaxies (see below). 

In this calculation, we included all the 15 fim members down 
to the limit of /(15^im) = 67^dy (i.e. 5.0a in the deep region and 
3.5a in the shallow region), but we found that this trend does not 
change if we limit the sample to relatively brighter 15 fim members 
(e.g. f(15fim) > 100/iJy which corresponds to 7.5a in the deep 
region and 5.2a in the shallow region). Also, as noted above, we in- 
cluded both resolved and unresolved 15 fim cluster members in this 
calculation. If we use only the resolved 15 ^m cluster members, the 
significance of the high fraction of optically red 15 ^m members in 
the medium-density environment becomes even higher. 

We now focus on the optically red 15 fim members (i.e. dusty 
red galaxies) that we have noted in Section |4~T1 We find 19 such 
galaxies in the entire field, and they are detected at significantly 
larger fraction in the medium-density environment compared to the 
low and high density environments (more than ~ 2 times larger). 
Detection in 15 fim of the optically red galaxies should mean that 
these galaxies still have on-going star formation activity, although 
their R — z' colours are as red as passively evolving galaxies. As 
described in Section |4~T| these galaxies are strong candidates for 
star forming galaxies heavily attenuated by dust. Our results there- 
fore suggest that the dusty star formation activity is induced by 
some mechanisms which are effective in the medium-density en- 
vironment. We plot the spatial distribution of these optically red 
15 cluster member galaxies in Fig. QT]on top °f a U the op- 
tically red phot-z members. It is quite noticeable that almost all 
dusty red galaxies are distributed exclusively along the filament 
including groups traced by the optical members, and in fact, this 
distribution is very similar to that of the member galaxies in the 
medium-density environment (see Fig. |7p). Although there are a 
few dusty red galaxies in the low and high density environments as 
well (Fig. I lot, almost all of them are located immediately outside 
of the filaments. 

To study the properties of the 15 members in the medium- 
density environment further, we investigate the z — L15 colours 
for the 15 fim members. The z'-band magnitudes approximate the 
stellar mass of galaxies (~ 5000 A in the rest frame for z ~ 0.8 
galaxies) and the 15 fim magnitude approximates the star forma- 
tion rate. Therefore, z — L15 colour approximately corresponds 
to "specific star formation rate (SSFR)" which is defined as star 
formation rate per unit stellar mass of galaxies. Since SSFR has 
a unit of inverse of time, it can be regarded as a timescale of star 
formation activity (i.e. galaxies which have high SSFR are consid- 
ered to have short star formation timescales). In Fig. [12] we plot the 
z — L15 colours for the resolved 15 fim members as functions of 
local density and z'-band magnitude. 

We notice that 9 out of 54 galaxies plotted in Fig. [12] have 
large colour indices of z — L15 > 3.7. We can consider that these 
galaxies are forming stars vigorously and efficiently for their stel- 
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Figure 9. Cumulative fractions of all the phot- ,2 members (solid line) and 
all the 15 fim members (dashed line) as a function of the distance from 
the cluster centre. The 15 fim members are much less concentrated in the 
central 1' region, while they dramatically increase within the ~ 1.5 — 2.0 
arcmin ring from the centre. These plots are made for galaxies with z' < 
22.5 only, because most of the 15 fim members are detected from such 
optically bright galaxies (see text in Section l4~TV 




Figure 10. Fractions of 15 fim member galaxies with f(l5fim) > 67/iJy 
as a function of the environment. The trend that the fraction is highest 
in the medium-density environment (especially for red galaxies) is to be 
noted. The open and hatched histogram in the upper panel represent the 
distribution of log S5 for all the phot-z members with z' < 22.5 and for 
15 fim members with z' < 22.5, respectively. The vertical dotted lines at 
log £5 = 1.65 and 2.15 represent the definitions of low, medium and high 
density environments. 
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Figure 12. z' — L15 colours as a function of the local density (a), and the z'-band magnitude (b). Filled and open symbols indicate optically red and blue 
resolved 15 fim cluster members, respectively. Galaxies above the horizontal solid line at z' — L15 = 3.7 are defined as "high SSFR galaxies". Note that 
these high SSFR galaxies are concentrated in the medium-density environment and concentrated at around 21.5 %, z' %, 22. The vertical dotted lines in 
panel(a) represent the definitions of the environment. The vertical dashed line and slanted solid line in panel(b) show the 5a limiting magnitudes in z' and 
L15, respectively. 
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Figure 11. The spatial distribution of the optically red 15 fim cluster mem- 
bers (open circles and triangles) on all optically red cluster members. The 
size of the open symbols indicates the 15 flux as in the previous figures 
(e.g. Fig. [8j. The large and small filled dots represent the z' < 22.5 and 
z' > 22.5 optically red cluster members, respectively. We can see that the 
distribution of optically red 15 fim members is aligned in the direction of 
the filament, and matches well with the distribution of the galaxies in the 
medium-density environment (groups/filaments). 



lar masses. This boundary roughly corresponds to 1 /SSFR ~ 0.25 
Gyr for z = 22 mag galaxies. This calculation is based on the 
stellar mass of a galaxy calcul ated from the z'-band ma gnitude us- 
ing the disk galaxy model in Koda mluBell. & Bowerl (19991) and 
the SFR calculated from the 15 fim magnitude. We hereafter call 
these galaxies "high SSFR galaxies". Interestingly, many of the 
high SSFR galaxies also tend to live in the medium-density en- 
vironment (see Fig. 112b). We note that 7 out of 23 (~ 30%) 15 fim 
members satisfy the high SSFR criterion in the medium-density 
environment, while only 2 out of 29 (~7%) galaxies in the low- 
density environments are high SSFR galaxies, and there is no such 
galaxy in the high-density environments. In Fig. l 12f b). we can see 
that the z'-band magnitudes of such high SSFR galaxies tend to 
be fainter than z' ~ 21.5. We show the spatial distribution of such 
high SSFR galaxies in Fig.[T3]with large open stars. The high SSFR 
galaxies also tend to be seen in or just around the filament, although 
their distribution may be slightly more extended than that of the 
dusty red star forming galaxies (i.e. open circles and triangles in 
Fig. O. 

We note that in Section I3.2l we found five possible AGN can- 
didates, at which three have X-ray detections and the other two have 
red N3— S7 colours. We find none of these five sources satisfy the 
high SSFR criterion of z' — L15 > 3.7. We find two optically red 
sources in these five AGN candidates (one in the low-density envi- 
ronment and the other in the medium-density environment). There- 
fore, we find that these AGN candidates do not strongly affect our 
results on the high SSFR galaxies and the dusty red galaxies. We 
stress here that even if these galaxies are really associated with pure 
AGNs and if we exclude these galaxies from our sample, our con- 
clusion is not changed. 

It is found that both the "dusty red galaxies" and the "high 
SSFR galaxies" prefer to live in the medium-density environment. 
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Figure 13. The distribution of the high SSFR galaxies (open stars). Small 
dots indicate all the phot-z member galaxies. The high SSFR galaxies 
are also distributed preferentially in the medium-density environment (i.e. 
groups/filaments). Open circles and triangles represent the "dusty red galax- 
ies", same as in Fig. 1111 

Since this environment is the one where colour distribution of 
galaxies start to change drastically from blue to red (see Sec- 
tion 14. 2\ , it is reasonable to think that these galaxies are related 
to such colour transition of galaxies, and they may be at the tran- 
sient phase from blue active galaxies to red passive galaxies during 
the course of hierarchical assembly to higher density regions. 



5 DISCUSSION 

In this section, we discuss our results in broader context of envi- 
ronmental dependence of galaxy evolution in clusters by compar- 
ing with previous works in the literature. Since our study is the 
first attempt to investigate around a z ~ 1 cluster with such a 
wide field coverage in MIR observation, we cannot directly com- 
pare our results with others for the region far out from the cluster 
centre. As for the cluster central regions, two z ~ 0.8 clusters 
(RXI0152 and MS 10 54, both at z = 0.83 ) ar e studied recentl y 
with Spitzer MIPS by IMarcillac et ail J2007t) and lBai et all d2007h . 
IMarcillac et ail d2007l) found that MIR cluster members are dis- 
tributed outside the two main cl umps traced by X-ray emissions in 
the RXI0152 cluster. Similarly. iBai et al.1 J2007h showed that there 
are few IR galaxies in the high-density regions of the cluster. For 
the RXI1716 cluster, we find in Section. 4.3 that the 15 /im mem- 
ber galaxies are distributed avoiding the central pa rt of the cluster, 
which is qu alitatively consis tent with the results of M arcillac et al.l 
d2007t) and IBai et al.l d2007l) . The absence of luminous infrared 
galaxies in the cluster centre is also reported for some lower- z clus- 
ters, inclu ding local Coma and Virgo clust ers by ISOCAM (see 
review by iMetcalfe. Fadda. & Biviancl [2005). Therefore, we have 
confirmed the low star formation activity in the high-density en- 
vironment such as cluster centres with the MIR observations to 
z ~ 0.8. It is natural to think that the red optical colours of the 



galaxies in cluster centres are primarily due to the lack of star for- 
mation activity and not due to the reddening by dust. 

We found that MIR bright cluster members are detected rela- 
tively far out from the clust er centre, especially in group/filament 
environment. Very recently, iFadda et al.l d2008h conducted a very 
wide-field MIR observation with Spitzer MIPS which covers two 
filaments around the Abell 1763 cluster (z — 0.23), and found 
that the fraction of starburst galaxies is more than twice larger 
in the filaments than in the inner region or outer fields of the 
clusters. The enhancement of star forming activity in filaments 
is also suggested in nearby clusters in the op tical studies (e.g. 
IPorter & Ravchaudhur"y||2007l ; |Porter et alj20"08r) . Our result for the 
cluster at z — 0.81 is qualitatively consistent with these studies 
in the sense that star forming activity is enhanced in the interme- 
diate density environment between low-density general field and 
the high-density cluster core. From the field study at z ~ 1, it 
is reported that the envir onment of LIRGs/UL I RGs is denser than 
that of field galaxies (e.g. lMarcillac et al.l2008l) . lEibaz et al.ld2007h 
studied a structure at z ~ 1 in the GOODS-North field with Spitzer 
MIPS, and showed that galaxies with strong star formation are pref- 
erentially seen in group centres. They suggested that, at z ~ 1, 
SFR of individual galaxies increases with increasing density up to 
a certain critical density and it decreases again a t higher density 
(see also ICooper et alj|2007h . For distant clusters, IPoggianti et al.l 
(2008) recently suggest a possible peak in the SFR-density relation 
based on the [Oil] line study of EDisCS clusters at z = 0.4 - 0.8. 
These studies are also qualitatively c onsistent wi t h our results, al- 
though the structures investigated in lElbaz et alj d2007h would be 
much poorer systems than the RXJ1716 cluster, judging from their 
weak X-ray detection. 

We now focus on the optically red star forming galaxies (i.e. 
"dusty red galaxies"). Dusty red star forming galaxi es in clus- 
ters w ere studied in some previous works. For example. ICoia et all 
d2005l) studied the CL0024+1654 cluster at z ~ 0.39 with ISO- 
CAM, and found that about half of the MIR sources in the cluster 
reside on the red-sequence. These dusty galaxies are not concen- 
trated in high-density regions. Base d on the COMBO- 17 data for 
the A9 01/902 clusters at z ~ 0.17, Wolf, Gray, & Meisen heimed 
d2005l) also showed that more than 1/3 of the red-sequence galaxies 
have dusty red SEDs. They also found that these dusty red galaxies 
prefer medium-density outskirts of the clusters and they are rare in 
low or high density environments. Since the redshifts of these clus- 
ters are lower than that of RXI1716 (z = 0.81), we cannot con- 
clude that these are the same populations as our dusty red galax- 
ies. However, judging from their optically red colours while hav- 
ing signs of on-going star formation activity, we expect that our 
dusty red galaxies are similar counterparts o f the dusty popu lations 
seen in the two low-z clusters stud ied by ICoia et al I d2005l) and 
IWolf, Gray. & Meisenheimerl d2005l) . As we showed in Fig. ll II our 
dusty red galaxies are preferentially located in the medium-density 
environment. It is interesting to note that these optically red star 
forming galaxies are common in the outskirt of clusters at all red- 
shifts through 0.1 S z g 0.8, although the number of cluster sam- 
ple is very limited. 

Finally, we discuss the physical mechanisms that can 
play major roles in galaxy truncation of star forming activity. 
IMarcillac etal] d2007l) found that MIR members in the RXJ0152 
cluster seem to be m ostly associated with t he in-falling late-type 
galaxies classified by Blakeslee et al. I d2006h . They suggest that a 
burst of star formation can occur during the galaxy in-fall process. 
Since we do not have spectroscopic data for our 15 /im cluster 
members, we cannot be sure at this stage that our 15 /im members 
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are physically associated to the in-falling galaxies, but it is likely 
that at least some of our 15 fim cluster members in the RXJ1716 
cluster are indeed the galaxies that are in-falling along the filament 
to the main body of the cluster by gravity. In fact, we have some 
good candidates of in-falling groups far away from the cluster cen- 
tre located at around (6'0, 0!0), (0!0, 6^0) and (~2'.0, 2'.0) in Fig. 
(8{a) as the iso-density contours clearly show. We show in Fig. ll4f a) 
an example of the candidates of an in-falling group located at (—0. 
'5, l'.5) in Fig.[8jb) which is ~ 1 Mpc away from the cluster centre. 
We need spectroscopic confirmation of the physical association of 
these systems to the cluster. 

We have revealed that the 15 fim members are found far out 
from the cluster core especially in the medium-density environment 
such as groups or filaments where the optical colour distribution 
strongly changes. Our results suggest that many galaxies which en- 
tered medium-density environment from the low-density field ex- 
perience starburst, and these galaxies are observed as dusty star 
forming galaxies and/or high SSFR galaxies around the RXJ1716 
cluster. Moreover, it would be natural to consider that the burst of 
star formation is linked to the colour and morphological changes of 
galaxies and truncation of star forming activity. The most likely 
physical mechanism at work in these medium density regions 
which also invol ves starburst would be galaxy-galaxy interaction 
or mergers (e.g. lHoplrinsetalJl2008h . We have two good exam- 
ples of 15 fim members which show prominent interacting features 
(Fig.ll4b.c1. The strong MIR emissions from these systems must 
be produced in the process of the galaxy interaction. These find- 
ings may suggest a link between the dusty star forming activity of 
galaxies and galaxy-galaxy interaction. However, it is difficult to 
firmly conclude that the majority of our 15 fim galaxies are acti- 
vated via galaxy-galaxy interactions or mergers. The seeing size of 
our ground-based optical image is ~0"7 and it is not sufficient to 
determine the morphology for all of our 15 fim galaxies and to find 
any interacting signatures. We should also admit that spectroscopic 
confirmation of membership of the phot-z members especially for 
the 15 /im members is crucial to draw any firm conclusions. Since 
our study reaches far out from the cluster core, it would have more 
contaminations than other studies limited in the cluster cores. The 
surface number density of the contaminant galaxies with z < 22.5 
mag is estimated to be ~ 0.93 arcmin -2 , using the galaxies in 
the control fields de fined in our Suprime-Cam field (see Fig. 3 of 
iKovama et ai]|2007h . A rough estimate shows that, among all the 
phot-2 members with z < 22.5 mag, ~10% and ~25% galaxies 
can be contaminant galaxies in the cluster region (i.e. high-density 
region) and the group region (i.e. medium -density region), respec- 
tively (see Fig. 4 of IKovama et al. I l2007l for the definitions of the 
cluster and the groups). Furthermore, we should keep in mind that 
the properties and/or fraction of the MI R galaxies may be di fferent 
from cluster to cluster, as suggested in iGeach et al.l ( |2006|) . even 
if we investigate clusters nearly at the same redshifts. Wide-field 
MIR observations on a larger sample of distant clusters is clearly 
required to obtain more general picture of the environmental effects 
on galaxy evolution. 



6 SUMMARY AND CONCLUSIONS 

We have performed a wide-field and multi-wavelength optical 
and infrared study of the distant galaxy cluster RXJ1716.4+6708 
(RXJ1716) at z — 0.81. A unique wide field coverage both in 
optical and infrared has enabled us to classify galaxies into three 
environmental bins, namely, high-density regions (cluster core), 
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Figure 14. Examples of the optical counterparts for 1 5 fim members (upper 
panels) and their 15 fim image (lower panels). Circles show the positions 
of 15 fim members and their radius are set to 8" in all the panels. Small 
squares represent the positions of phot- ,2 selected member galaxies and the 
size of the squares are 4" X 4" . (a) is a candidate for an in-falling group, 
(b) and (c) are the prominent cases of the interacting systems with 15 fim 
emissions. 



medium-density regions (cluster outskirts, groups, filament), and 
low-density regions (field), and has thus allowed us to investigate 
galaxy properties as a function of environment along the structures 
in and around the distant cluster. 

We find many of the 15 fim cluster members show interme- 
diate optical colours between the red sequence and the blue cloud. 
This may indicate that these 15 fim cluster members are actively 
star-forming galaxies but attenuated by dust hence showing in- 
termediate optical colours, although relatively fewer detection at 
15 fim of blue galaxies may be partly because they tend to be opti- 
cally faint. 

We quantified the environment around the cluster using the 
local projected number density of cluster member galaxies, and 
confirmed that the optical colour distribution starts to dramatically 
change at the "medium" density environment that corresponds to 
groups and/or filaments. We showed that the 15 fim members are 
very rare in the high-density cluster centre. This is probably due 
to the low star forming activity in such regions. However, inter- 
estingly, the fraction of the 15 /mi-detected cluster members in the 
medium-density environment is as high as in the low-density fields, 
despite the fact that optical colours of galaxies start to dramatically 
change from blue to red in the medium-density environment. Al- 
though the statistic is not very good, the fraction is slightly higher 
in the medium-density environments even compared with the low- 
density fields. 

We also find that dusty red galaxies (optically red 15 /mi 
cluster members) and the galaxies with high specific star forma- 
tion rates (red in z' — 15 /mi colour) are both concentrated in the 
medium-density environment. These results may suggest that the 
star formation activity in galaxies is once enhanced by some phys- 
ical processes which are effective in group/filament environment 
(e.g. galaxy-galaxy interaction), before their star forming activity 
is eventually truncated and they move on to the red sequence. 

We stress that all these new findings are based on the widest 
field MIR observation of a z ~ 0.8 cluster so far. There is no other 
study that covers such a wide field around clusters in MIR at z > 
0.8. Since our study is a case study for just one cluster at z = 0.81, 
we are desperately in need for a larger sample of distant clusters 
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viewed in the infrared regime and at the same time covering a wide 
field of view so that we can witness the galaxy truncation in action 
in the in-fall regions of distant clusters along the filaments. This is 
essential in order to confirm this interesting trend and to obtain a 
general view of galaxy evolution. 
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